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We apply the concepts of nonlinear guided-wave optics to a Bose-Einstein condensate (BEC) 
trapped in an external potential. As an example, we consider a parabolic double-well potential 
and derive coupled-mode equations for the complex amplitudes of the BEC macroscopic collective 
modes. Our equations describe different regimes of the condensate dynamics, including the nonlinear 
Josephson effect for any separation between the wells. We demonstrate macroscopic self-trapping 
for both repulsive and attractive interactions, and confirm our results by numerical simulations. 
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A system of interacting bosons confined within an ex- 
ternal potential at zero temperature can be described by 
a macroscopic wave function having the meaning of an 
order parameter and satisfying the Gross-Pitaevskii (GP) 
equation W]. The GP equation is a nonlinear equation 
that takes into account the effects of the particle inter- 
actions through an effective mean field, and it describes 
the condensate dynamics in a confined geometry. Similar 
models of the confined dynamics of macroscopic systems 
appear in other fields, e.g. in the case of an electron 
gas confined in a quantum well, or optical modes of a 
photonic microcavity |2J. In all such systems, confined 
single-particle states are restricted to discrete energies 
that form a set of eigenmodes. 

The physical picture of eigenmodes remains valid in the 
nonlinear case |^, and nonlinear collective modes corre- 
spond to the ground and higher-order (excited) states of 
the Bose-Einstein condensate (BEC) |^. Moreover, it 
is possible to observe at least the first excited (antisym- 
metric) collective mode experimentally H, through the 
collapses and revivals in the dynamics of strongly cou- 
pled two-component BECs Iq]. The interest in the non- 
ground-state collective modes of BECs has grown dra- 
matically with the study of vortex states, very recently 
successfully created in the experiment . 

The modal structure of the condensate macroscopic 
(ground and excited) states allows us to draw a deep 
analogy between BEC in a trap and guided-wave optics, 
where the concept of nonlinear guided modes is widely 
used [p[. The physical description of confined conden- 
sate dynamics in time is akin to that of stationary beam 
propagation along a nonlinear optical waveguide, with 
the BEC chemical potential playing the role of the beam 
propagation constant. As is well known from nonlinear 
optics, the guided waves become coupled in the presence 
of nonlinearity, and the mode coupling can lead to the 
nonlinear phase shifting between the modes, power ex- 
change, and self-trapping. 

In this paper, we apply the concepts of nonlinear 
guided-wave optics to the analysis of mode coupling and 
intermodal population exchange in trapped BECs. As the 
most impressive (and also physically relevant) example 
of the applications of our theory, we consider the BEC 



dynamics in a harmonic double-well potential, recently 
discussed in the literature |^ . We study the coupling be- 
tween the BEC ground-state mode and the first excited 
(antisymmetric) mode in such a potential, and derive the 
dynamical equations for the the complex mode ampli- 
tudes, valid for any value of the well separation. Our 
model comprises, in the limiting case of large separation, 
the theory of Josephson tunneling developed for weakly 
interacting condensates in two separate harmonic traps 
|lO| . In the limit of close separation, our theory describes 
a nonlinear population exchange between the interact- 
ing modes, similar to the effective Rabi oscillations in 
two-component BECs, studied both theoretically M and 
experimentally |ll| ]. 

We consider the macroscopic dynamics of BEC in a 
strongly anisotropic external potential, U — ^muj^{Y^ + 
Z^ + XX^), created by a magnetic trap with a character- 
istic frequency to. In the case of the cigar-shaped trap, 
A ^ 1, the collective BEC dynamics can be described by 
a one-dimensional GP equation. Details of the derivation 
and normalization can be found, e.g., in Refs. |12[| . The 
GP equation for the longitudinal profile of the normalized 
condensate wave function takes the form: 
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C/(x)V' + cr|V'lV = 0. 



(1) 



According to normalization, the number of the conden- 
sate particles M is defined as A/" = {Tiuj/2Uo^/\)N , where 
Uq = 47r?i (a/m) characterizes two-particle interaction 
proportional to the s-wave scattering length a, and the 
functional N = J_ \ip\'^dx is the integral of motion 
for the normalized nonstationary GP equation (w). The 
value of CT = — sgn(a) = ±1 in front of the nonlinear term 
is defined by the sign of the scattering length of the two- 
body interaction, repulsive for a > 0, and attractive for 
a < 0. The potential U{x) = k{\x\ — xq)^ (below we take 
fc = 1) describes the double- well structure of the trap in 
the longitudinal direction. 

In the linear limit (i.e. for an ideal non-interacting 
gas), we should consider cr = 0, and the exact station- 
ary solutions of Eq. (||) in the form t/^ix, t) — ^j{x)e^^^* 
are found in terms of the parabolic elliptic functions |l3] 
that define a set of confined stationary states existing 



at certain discrete values of Pj (linear eigenmodes). For 
(7 ^ 0, we can introduce, in a similar way, a set of non- 
linear eigenmodes Q described by real functions ^j{x) 
that satisfy the following nonlinear equation. 



dx^ 



P,^,-U{x)^, 
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0. 
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While in the case ct = the eigenvalue for each mode f3j is 
unique for any given trap separation xq , in the nonlinear 
case there exist families of localized solutions $j charac- 
terized by the dependence of the norm Nj = J dx^^ (x) 
on /3j, and the eigenvalue now becomes a parameter of 
a continuous family f^. Figure 1 shows two examples 
of the ground-state mode ^o{x) and first-order excited 
mode $i(a;) of the BEC with A^o = ^i but /3o 7^ /3i, for 
different values of xq. The dependencies /3o and /3i on 
the trap separation xq are quite different for two signs of 
CT, as is seen in Fig. 2. 




FIG. 1. Confining potential (bold) with the ground (solid) 
and first excited (dashed) collective modes for (a) xq — 0.6, 
and (b) xo = 2.5 (a = -1, A^o = A^i = 5.0). Dotted 
lines - corresponding values of the chemical potential: (a) 
/3o = 1.974, /3i = 3.162, and (b) f3o = 1.889, /3i = 1.925. 




FIG. 2. Dependencies /3o(a;o) (lower curve) and /9i(xo) 
(upper curve) for cr = (dotted) and a — ±1 (solid, 
No = Ni = 5.0). Dashed - results of the variational approach. 



where bj [j = 0, 1) are the complex amplitudes, and 
$j(a;) may, in general, be any two solutions of Eq. (||). 
Then, the BEC dynamics can be deduced from the rate 
equations for the modal amplitudes bj. A similar ap- 
proach has previously been employed in Refs. pO| | to 
describe coherent tunneling between two largely sepa- 
rated harmonic traps, with the basis $j consisting of 
the ground-state modes of individual potential wells. In 
contrary, our basis eigenfunctions are the local nonlinear 
modes of the entire double-well trap which can be found 
exactly for any given trap separation. Our approach is 
therefore similar to the analysis of the power transfer 
between the cores of a nonlinear optical coupler usually 
carried out in the terms of the supermodes, i.e. the local 
modes of a composite core Q . 

To derive the equations for the complex mode ampli- 
tudes &j(i), we use the standard procedure, substituting 
the ansatz (0) into the nonstationary GP equation (1^), 
and averaging over the spatial dimension after multiply- 
ing the GP equation by either <^o{x) or <I'i(x). This yields 
a system of two coupled equations. 
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and the 
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/dx<i>f(x)$2(:,) 



Here Vt = 2(/3i - /3o) - 2a{CaNo ~ 
coupling coefficients are defined as 

Coi = 7oi/(^o^^i), where 7ij = 
and Bj are the normalized mode amplitudes, Bj(t) 
y/N~jbj{t) exp{—iaCjNjt). These equations conserve the 
total norm |Bop-|-|i?ip = no{t)-\-ni{t) = n, where no and 
ni have the meaning of the time-dependent population 
numbers for the two macroscopic states, and n = N = 
^ol^oP + .^il^iP- It is important to note, that the form 
of the rate equations does not depend on the normaliza- 
tion conditions for the basis functions. For example, the 
condition /_ \^j\^dx — J_ \ip\'^dx — I simply imposes 
the constraint \Bj\ = \bj\, so that n^ N ^ |6op + |5ip. 
Separating the amplitudes and phases as Bj{t) — 
yjnj{t) exp[— i^j(i)], we obtain a system of coupled equa- 
tions for the population difference of the two states A — 
iT-i — no, and the relative phase shift, Q = 2((/>o — (t>i) — ^t, 



a{n^ -A^) sine. 



dA 

~dt 

^ = -5 + CT(Co + Ci)A-2crCoi(2 + cose)A, 
dt 



(4) 



To develop a coupled-mode theory for BECs, we con- 
sider the mode interaction in a double- well potential and 
assume that the condensate wave function is described 
by a superposition of two modes of different symmetry, 
i.e. symmetric and antisymmetric, 

V'(x,t) =6oW<i>o(2:)e-^''°* + 6i(t)<i>i(x)e-'^i*, (3) 



where 5 — 2(/?i — 
System (||) can be 
-dH/de, de/dt = 
(T(n^- A^)Coi cos 
mechanical analog 
tion of a non-rigid 
and a generalized 



Po) + cj[in - 2iVo)Co - (n - 1Ni)Cx\. 
rewritten in a canonical form, dA/dt = 
= dH/dA, with the Hamiltonian: H ~ 
e + a [{Co + Ci)/2 - 2Coi] A'^-SA. A 
;y of this system may describe the mo- 
pendulum with angular momentum A 
angular coordinate 8. On the other 



hand, Eqs. (^ closely resemble the dynamic equations 
for the guided power and phase difference of two nonlin- 
early interacting orthogonally polarized optical modes in 
a birefringent fiber [g . The exact solution of the system 
(0) can be obtained in terms of elliptic Jacobi functions 
and will be presented elsewhere. 

The dynamics described by Eqs. (0) depends crucially 
on the values of the coupling coefficients Co, Ci, and 
Coi , which are determined by integration over the eigen- 
mode profiles, so that the results can be quite different 
for the two signs of a. Moreover, the condensate dynam- 
ics changes with the separation of the potential wells, 
as governed by the dependencies C[j{xq), Ci{xq), and 
Coi(a;o), which differ drastically for cr = ±1 (see Fig. 
3). As expected from the linear theory fl^ and results 
for cr = — 1 PI, the energy spectrum becomes degener- 
ate for large separation (see also Fig. 2), and the cou- 
pling between the collective modes becomes more coher- 
ent. Importantly, for a — +1, this happens at the values 
of separation smaller than those for cr = — 1. 






FIG. 3. Coupling coefficients as functions of separation for 
(a) a = -1-1, and (b) a = —1. Dashed - results of the varia- 
tional approach. 

The routine of calculating the coupling coefficients nu- 
merically can be bypassed by employing the variational 
approach^ using the trial functions in the form of a lin- 
ear superposition of the ground states of isolated traps: 
*o,i = ^o,i{exp[(2: - xo) V2ao,i] ± exp[(x -f xq) V2ao,i]}' 
and the Lagrangian of the stationary GP equation. 
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By inserting the trial functions into the corresponding 
variational integral, an explicit, although algebraically 
complicated integrated Lagrangian is obtained. The vari- 
ational equations with respect to the parameters v4o,i 
and ao.i yield relations which determine the characteris- 
tic eigenvalues Pq^i{xq) and the coefficients Co, Ci, and 
Coi . These relations can be further simplified in the limit 
xo ^ 1, but must be solved numerically for a general 
case. Comparisons between the variational predictions 
and the results obtained by solving Eq. (g) numerically 
are shown in Fig. 2, for /?o,i(a;o), and in Fig. 3, for Co, 
Ci , and Coi . The agreement is seen to be satisfactory. 



To visualize the population dynamics, in Fig. 4 we plot 
the phase portraits {O, A} of the dynamical system (Q) 
for the case a — +1 and Nq — Ni. For convenience, the 
population difference. A, is measured in the units of n. 
For small separations [Figs. 4(a,b)], while the coupling 
coefficients are sufficiently different, there are only two 
fixed states of the relative population: A = ±n, which 
corresponds to either no = or ni = 0. In both these 
states, the phase is unbounded, i.e. it is a linear func- 
tion of time. The other phase trajectories in Figs. 4(a,b) 
represent the dynamical states with the running phase, 
that is a delocalizcd phase. The mechanical analogy of 
this phenomenon is simple ||I3 : it corresponds to a self- 
sustained steady closed-loop rotation of a non-rigid pen- 
dulum around its support. In terms of the condensate 
dynamics, these states describe the nonlinear Rabi-type 
oscillations between the ground and first excited macro- 
scopic states, for small xq, and the Josephson-type tun- 
neling between the two potential wells, for a sufficiently 
large separation. Remarkably, the population of either 
well is never completely depleted. 
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FIG. 4. (a)-(d) Phase plane A(e) for Eqs. (|) at cr = +1, 
N = n = 1, and xq = 0.0, 1.5, 1.7, and 3.0, respectively. 

The phase plane shown in Fig. 4 also reveals the 
existence of macroscopic quantum self-trapped (MQST) 
states, predicted and described in llQ] for weakly interact- 
ing BECs in largely separated traps. The MQST states 
are characterized by a nonzero average population imbal- 
ance. As the separation between the wells grows, a bifur- 
cation of the fixed points on the phase plane {A, Q} oc- 
curs, and the stable centers, corresponding to the MQST 
states with a trapped phase, appear [see Figs. 4(c,d)]. 
This occurs at a certain a;g'', for which the condition 
S — n(6Coi — Co — Ci) is satisfied. For cr = +1 and 
n = 1, for example, Xq^ ~ 1.48 which agrees very well 
with the corresponding result of the variational approach, 
Xq^ ~ 1.41. For larger separation, when Cq ~ Ci '--^ Coi, 
the positions of the centers are approximately given by: 
A « (/?o — /3i)/(2Coicr) at = ±(2m)7r, where m is inte- 
ger. With increasing separation, as (/3o — (3i) ^ 0, these 



fixed centers move towards the line A = 0, and the sad- 
dles form between them, so that the MQST states, other 
than those identical to the ground states of the individual 
wells, cease to exist [see Fig. 4(d)]. The oscillations of the 
population imbalance around the stable fixed points with 
the trapped phase have been identified in ||l^ as 7r-states. 
Clearly, the effect we observe here is qualitatively similar, 
except for the average value of the trapped phase which, 
due to a different choice of the basis eigenfunctions and 
definition of Q, is equal to zL{2m)Tr. 
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FIG. 5. (a) Phase trajectories A(0) calculated from Eqs. 
for a = —1, xq = 1.7 (thin) compared with the numerical 
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(I) 

solution of the GP equation (thick), (b) \tp\^ in the MQST 
in (a), at the normalized times t = (solid) and 



t = 100 (dashed). 



i/)| corresponding to the point 2 in 



(a) at i = (solid), and corresponding to the point "3" in 



at t = 20 (dashed), 
by a dotted curve. 



Double-well potential is shown in (b,c) 



To compare the predictions of our coupled-mode the- 
ory with the actual dynamics of the BEC in a double- 
well potential modeled by the GP equation, we solve 
Eq. (H) numerically. As an initial condition, for both 
a = ±1, we take ip{x,0) = bo{0)^o{x) + bi{0)^i{x), with 
^o(0) + &?(0) = l/No = 1/A^i. In Fig. 5(a), the phase tra- 
jectories A (9) are compared with those calculated using 
Eqs. (0), for a — —1, and the trap separation corre- 
sponding to sufficiently dissimilar values of (ij and Cj 
(see Figs. 2 and 3). It is clear that the approximate 
equations of the coupled-mode theory correctly describe 
the dynamics of the condensate in the states with a run- 
ning phase [see Fig. 5(c)], as well as the position of the 
MQST states, one of which is shown in Fig. 5(b). Per- 
forming a similar comparison for different a and xq, we 
can conclude that the eigenfunctions ^j{x) represent a 
good basis for the modal decomposition of the macro- 
scopic condensate wave function tp{x,t). The adiabatic 
evolution of the eigenfunctions with time, although lead- 
ing to slight deviations of the condensate states from the 
exact MQSTs [see Fig. 5(a)], does not introduce signif- 



icant damping into the system, and therefore does not 
lead to a dramatic switching between the states. 

In conclusion, we have employed the concepts of the 
nonlinear guided-wave optics and developed, for the first 
time to our knowledge, a consistent coupled-mode the- 
ory for BECs. We have studied the BEC dynamics in 
a double-well harmonic trap, and verified the results by 
numerical simulations of the nonstationary GP equation. 
The strong advantage of our theory is its ability to de- 
scribe the condensate dynamics for any well separation, 
including the Josephson tunneling effect at large sepa- 
rations, mode coupling and Rabi oscillations in a single 
harmonic well, and the macroscopic self-trapped states 
in the crossover regime. 
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